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Abstract 

In this work, we calculate the next-to-Ieading-order (NLO) QCD corrections to the process 
p + p — > J /ip + ^ via the color-singlet mechanism at the LHC. The results show that the partial 
cross section {p'l^^ >3GeV and [y^/'^'Tj <3) is enhanced by a factor of about 2.0, and the differential 
cross section can be enhanced by two orders of magnitude in the large transverse momentum region 
of J /ip. Furthermore, the polarization of J/ip changes from transverse polarization at leading-order 
to longitudinal polarization at NLO. For the inclusive J /if: hadroproduction, it is known that the 
color-octet contributions are one order of magnitude larger than the color-singlet contribution, 
and the polarization distribution is dominated by the color-octet behavior at NLO. In contrast, for 
J/V'+T production the color-singlet contribution is of the same order as the color-octet contribution, 
and the polarization distribution arises from both the color-singlet and color-octet. Therefore, 
measurements of J /ip production associated with a direct photon at the hadron collider could be 
an important supplement to testify the theoretical framework treating with the heavy quarkonium. 
In addition, an NLO QCD correction to T -|- 7 production at the LHC is also presented in this 
paper. 

PACS numbers: 12.38.Bx, 13.25.Gv, 13.60.Le 
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For a long time, the productions and decays of heavy quarkonium have been an ideal 
laboratory to investigate quantum chromodynamics. The large mass of the heavy quarks 
provides a large scale for production and decay processes. This large scale makes the factor- 
ization of the calculations for these processes available, and the hard part in the calculations 
can be calculated perturbatively. Conventionally, the color-singlet mechanism(CSM) [1] was 
used to describe the production and decay of heavy quarkonium. However, the CSM has 
encountered many difficulties in various theoretical 2I and experimenta^ aspects [sl. In 1995, 
non-relativistic quantum chromodynamics(NRQCD) was put forward[4]. By including the 



contributions from the high Fock states, NRQCD can overcome the theoretical difficulties 



iesfl. 



At the same time, the NRQCD predictions are consistent with experimental data[6|. Al- 
though NRQCD has had many successes, there are also many problems in the production 
of heavy quarkonium[7t]. 

In 1995, Kramer calculated the NLO QCD corrections to inclusive J /if) photoproduction 
and found that the data from HERA could be understood by just including the high order 
QCD corrections in CSM 8]. Recently, many studies [sl, [lO] have shown that NLO QCD 
corrections in CSM to J/ ip related production processes at B-factory change the theoretical 
predictions dramatically. Furthermore, at the Tevatron, the NLO QCD corrections in CSM 



enhance the total cross section by a factor of 2 and give large values for the transverse 
momentum pt distribution of J/ip^ especially in the large pt region |ll. [l^. The polarization 
status of J/ip drastically changes from transverse-polarization dominant at leading-order 



(LO) into longitudinal-polarization dominant at NLO for J/ip hadroproduction 12 



laj. Ar- 



toisenet and his collaborators suggested that the pt distribution of T(IS') at the Tevatron 
can be interpreted in the CSM by taking the NNLO real part into account [14] . It is quite 
clear that higher order QCD corrections play a very important role in theoretical predictions 
for the production of heavy quarkonium, and their effects should be carefully investigated 
in all the relevant processes. 

In the framework of CSM, the associated production of J/^/^ + 7 at a hadron collider 
was first proposed as a good channel to investigate the gluon distribution in the proton 
with a relatively clean signal[l^. Soon thereafter, it was used to study the polarized gluon 
distribution[16] and the production mechanism of heavy quarkonium 17]. In reference [l8|. 
the associated production of J/?/) + 7 at the Tevatron has been considered in the CSM at 
the LO, and the results show that the contribution from gluon fusion process is dominant 
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over that from the fragmentation processes. Kim investigated the contribution of color-octet 
processes to the hadroproduction of J/^/'-|-7 and found that the color-octet contributions are 



dommant in the la^ge p, reg,o„M. The eontr.but.ou fro,„ fragmentation p_ is smaller 

;han the fusion contribution within the NRQCD framework at the LHC[20|j. In reference 



2l| . the authors obtained a theoretical prediction for J/-?/; + 7 hadroproduction at the 
hadron collider, and the numerical results show that the transverse momentum distribution 
of J /if) production is smaller in the CSM than that in the color-octet mechanism (COM) 
in the large pt region at LO. To further study the effect of the NLO QCD corrections on 
heavy quarkonium hadroproduction, in this paper we calculate the NLO QCD corrections 
to J/^ + 7 hadroproduction at the LHC and present theoretical predictions for the pt 
distribution of the production and polarization for J /if). 

We employ the automated Feynman Diagram Calculation package (FDC) to perform 
the analytic evaluation of all the processes. FDC is a powerful program based on the 



LISP 



anguage designed to automate the NLO calculation and was initially developed by 



Wang|22l|. Its one-loop part was recently completed by Wang and Gong 
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It has recently 
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24 



25|. 



been successfully applied to several quarkonium production processes 

For the process p + p^J/ip + ^aX LO, only the gluon fusion process g + g ^ J/ip + 7 
contributes at the partonic level with six Feynman diagrams, which is similar to that of 
g + g—>-J/ilj + g in inclusive J/ip hadroproduction. The LO total cross section is obtained by 
convoluting the partonic cross section with the parton distribution function (PDF) Gg{x, fif) 
in the proton: 

cr-^ = y dxidx2Gg{xi,fXf)Gg{x2,fXf)a^, (1) 

where /xj is the factorization scale. In the following, a represents the corresponding partonic 
cross section. 

When the cross section is expanded to NLO on as, the coupling constant of quantum 
chromodynamics, there are one virtual correction and three real correction processes. They 
are listed as the following: 

g + g^J/^ + -f, (2) 
g + g ^ J/^ + j + g, (3) 

qiq) + g ^ Jf-ip + i + qiq), (4) 

q + q-^ J/^ + ^ + g. (5) 



In calculating the virtual corrections in Eq.([2]), there are 111 Feynman diagrams at NLO. The 
ultraviolet (UV) and infrared (IR) singularities are normalized and separated by using the 
dimensional regularization. Furthermore, the same renormalization scheme as in Ref. [l2.[l3| 
is applied to redefine the quark mass, coupling constant and the quark or gluon fields, 
and then the renormalized amplitude without UV singularities is obtained. The Coulomb 
singularity, which comes from the diagrams with a virtual gluon connecting the quark and 
anti-quark pair in J /if), is regulated by introducing a small relative velocity v between the 
quark and anti-quark pair and is absorbed into the redefinition of the wave function of J /ip. 
Therefore, the one-loop amplitude is free of UV and Coulomb singularities and the 
virtual corrections to the NLO cross section are expressed as 

^ oc 2Re(M^M^*), (6) 

where is the amplitude at LO. It is UV and Coulomb finite, but has IR singularities. 
To obtain an infrared-safe cross section, it is needed to cancel the IR singularities by adding 
the contributions from real processes at NLO. 

For the real processes in Eqs.Q, (jl]) and ([3]), there are IR singularities in the phase space 
integration which are divided into soft and collinear singularities. The soft singularity is from 
a soft gluon emitting from the initial gluons. The collinear singularity is from the initial 
gluon (or quark) emitting a gluon which is nearly parallel to the parent particle. It is easy to 
find that soft singularities caused by emitting soft gluons from the charm quark- ant iquark 



pair in J / are canceled by each other. Using the standard two cut-off slicing method 26|] , 
we decompose the phase space into three regions by introducing two small cutoffs 6s and 6c- 
Therefore, convoluting the parton level cross section with the parton distribution function, 
the cross section of the real processes is represented as 

a^ = a^c^^s^^HC^^HC_ (7) 

Here cr^*" is the contribution from the hard noncoUinear part of the phase space that is 
IR finite, is the soft part and a^^ is the hard collinear part. After absorbing the mass 
factorization parts into the PDFs, there is a finite term remaining which is similar to 
that in Eq.(43) in Ref. [131]. By adding all the contributions together, a finite total cross 
section is obtained. 

The polarization of J/ip is an important physical measurement to identify different pro- 
duction mechanisms. Therefore, the polarization is calculated at NLO. Here we describe the 



10 




.. I .... I .... I .... I .... I .... I .... I .... I .. . 

5 10 15 20 25 30 35 40 45 50 
P, of .i/A- (GeV) 



FIG. 1: Transverse momentum distribution of J /if) production with = fJ-f = V (2mc)^ +Pf at 
the LHC. 

decay of J/iJj in its rest frame with the 3- momentum direction of J/ip as the z-axis. The 
polarization parameter a is defined as 

dar/dpt - 2daL/dpt 



(8) 



dar/dpt + 2daL/dpt ' 

where cr^ and are the cross sections of transversely and longitudinally polarized J/ijj 
respectively, a — —1 corresponds to fully longitudinal polarization and a = 1 to fully 
transverse polarization. 

To check the gauge invariance for each subprocess, the polarization vectors are replaced 
by the corresponding momentum in the numerical calculation and the gauge invariance is 
obviously observed at double precision level. Since the two phase space cutoffs are chosen 
to handle the IR singularities of the real processes, we varied the cutoffs by a few orders 
of magnitude to numerically check the independence of the results from the cutoffs and 
obtained consistent results within the error tolerance. 

In the numerical calculation, we use q;=1/137, rric—l-BGeV, and Mj/^ — 2mc. Prom the 
leptonic width of J/V', the wave function at the origin is extracted by 



Ml, 



(9) 
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We obtain \R^^{0)\'^=OMAGeV^ at a, = a,(Mj/^) = 0.26, T^e = 5.55keV. The Cteq6Ll 
and Cteq6M[27!] are used in the calculations at LO and NLO respectively, with the corre- 
sponding as running formula in Cteq6 being used. The renormalization scale fir and the 
factorization scale /ij are set to yUj. = /i/ = a/ (2mc)^ + p^. The center-of-mass energy at 
the LHC is 14TeV and the rapidity cuts for both J/^/'(T) and the direct photon are chosen 
as \y\ < 3. Because the perturbative expansion does not work well in the small transverse 
momentum region and the large rapidity region, we impose pt and y cuts on the numerical 
calculation. In these processes, the direct photon has nothing to do with the infrared sin- 
gular structure, so it can be identified with any experimentally acceptable kinematic cuts. 
Our numerical calculation is obviously consistent with this expectation. Of course, the cut 
conditions imposed on the direct photon will affect the numerical results. In the following, 
we fix the rapidity cut of the photon as |y| < 3 and vary its pt cuts to calculate the numerical 
results. 

By replacing the corresponding parameters as 

rric ^ rUb, Mj/^ ^ Mr, Cc ^ Cb 
R^{Oy/^ ^ R,{0)^ ^ n/ = 3^n/ = 4, (10) 

the results of T can be obtained. For the bottom quark mass and the wave function of T at 
the origin, mfe=4.75GeV and |i?T(0)P=7.74GeV^ are used. We use the Eq.Q to calculate 
the \Rr{0)\^ with ^.(Mt) = 0.18 and Tx^ee = 1.34keV. 

The Pt distribution of J/ip with different pt cut conditions on the associated photon is 
shown in Fig. [H It is noteworthy that there are two pt cut conditions, of which one is the 
J/ ip Pt cut condition and the other is the associated photon pt cut condition since both J/ ijj 
and the associated photon have to be measured experimentally. The NLO numerical results 
of the Pt distribution of J / ip are larger than the LO results by 1 ~ 2 orders of magnitude in 
the large pt region, and different cut conditions only affect the results in the region near the 



minimum endpoint in pt. In the case of inclusive J/ip hadropro duct ion at NLO in COM 24| . 
the results show that the NLO QCD corrections change the results by about 10 percent. 
Therefore, we can expect that the NLO QCD corrections will modify the results of J/ip + 'y 
at NLO in COM slightly, since the process is quite similar to inclusive J/ ip hadroproduction 
at NLO in COM. Fig. [2] shows that the NLO QCD corrections also change the pt distribution 
of T, but with less enhancement than that of J/ip. 
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FIG. 2: Transverse momentum distribution of T production with fir = fJ-f = \J (2rri;,)^ + at the 
LHC. 

The theoretical predictions of the J j^l) polarization at NLO are presented in Fig. [3] with 
different cut conditions. They are similar to the result of J/^^ inclusive hadroproduction. 
The J I'll) polarization drastically changes from the more transverse polarization at LO to a 
more longitudinal polarization at NLO. Also, the polarization parameter a oi J/ip at NLO 
becomes closer to -0.9 as pt increases. In the figures, the different curves are plotted with 
different pt cut conditions. It is clear that the pt cut condition for the associated photon 
has a greater effect near the minimum endpoint of J/ip pt and a smaller effect in the large 
Pt region. From Fig. HI the polarization of T also changes from transverse-polarization 
dominant to longitudinal-polarization dominant, and a at NLO becomes closer to -0.8 as 
Pt increase. In the small pt region, the convergence of the numerical calculation becomes 
worse, and it requires too much CPU time to improve the precision of the calculations. At 
the same time, the denominator of a changes its sign in the neighborhood of a pt point in the 
small Pt region, meaning that the value of a undergoes a vast change around that point. The 
small value of the denominator also amplifies the error in the calculation, indicating that 
the convergence of the perturbative expansion is bad and that higher order contributions 
are important in the small pt region, especially for the pt distribution of the polarization 
parameter a. Therefore, the small pt region is further discarded in our presentation of a. 
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FIG. 3: Transverse momentum dependence of J/V' polarization with fir = /J-f = \/ (2mc)^ + at 
the LHC. 

In Fig. [5], the /i dependence of the partial cross section {p'l^^^'^^ >3GeV,|?/| < 3) for 
J/iIj{T) +7 + X production at LO and NLO are shown with our default choice, fir = fJ'f = fJ'- 
In other words, the factorization scale equals to the renormalization scale, and fi ranges from 
/io/2 to 8yUo- There is no pt cut for the associated photon. It is obvious that QCD corrections 
reduce the fi dependence of the partial cross section for J/ip production and moderate the 
fi dependence for T production. 

To illustrate the influence of the cut conditions on the partial cross section, we present 
LO and NLO partial cross sections for J/^/'(T) +7 + X production at the LHC with different 
cut conditions in Table [H In Fig. El we plot the dependence of the partial cross section on 
the Pt cut of the direct photon with the pt cut of J/ipiT) fixed at lOGeV. The plots show 
that the contributions from small pt photon are large and the partial cross section decreases 
as the Pt cut on the photon increases. For the production of T, when the pt cut on the 
photon becomes too close to the pt cut of T, the partial cross section decreases rapidly and 
the influence of the error becomes severe. In this case, too many contributions from the real 
processes are cut off and the perturbative calculation becomes worse. Therefore, the partial 
cross section cannot be taken seriously when the pt cut condition of the photon approaches 
the Pt cut condition of T. 
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FIG. 4: Transverse momentum dependence of T polarization with ji^ = Hf = \J {2mb)'^ + at the 
LHC. 

In summary, we have calculated NLO QCD corrections to the production of J/iIj{T) as- 
sociated with a direct photon at the LHC in CSM. For J/ip production, the partial cross 
section is enhanced by a factor of about 2.0 with J/ip transverse momentum cut pt >3GeV 
and rapidity cut \y\ < 3 for both J/ip and the direct photon. The transverse momentum 
distribution of J/tp at NLO is enhanced by 1 ~ 2 orders of magnitude over LO calcula- 
tions as Pt becomes larger. The J/ip polarization is calculated, and the results show that 
the J/tp polarization drastically changes from transverse-polarization dominant at LO to 
longitudinal-polarization dominant at NLO. The situation is quite similar to the case of 
NLO QCD corrections to the inclusive J/ip hadroproduction. It can be seen that NLO 
results on J/iJj + '-f hadroproduction in CSM is of the same order of magnitude as LO results 



in COM 



19 



2l[ |. As a reasonable estimate from the experience of NLO QCD correction cal 



culation to inclusive J/tp hadroproduction in COM 2J], the NLO pt distribution for J/ip + 'j 
hadroproduction in COM will be changed slightly and the polarization will remain almost 
the same. Therefore, the results of J/ip production associated with a direct photon from the 
CSM and COM are of the same order of magnitude at NLO. Therefore, the J/tp polariza- 
tion measurement would be able to distinguish the contributions from the CSM and COM 
since the theoretical predictions of J/ ip polarization are obviously different for the CSM and 
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FIG. 5: The dependence of the partial cross section (p/''^^^'* >3GeV) with iJ,o = 

COM. In contrast to the case of J/'^ inclusive hadroproduction, the NLO result of CSM is 
about one order of magnitude smaller than the NLO result of COM, and so there is no way 
for the CSM part to play an important role in the final J/ip polarization distribution. 

To measure the production of J/iIj{T) associated with a direct photon at the LHC, the 
direct photon must be identified. To cut the enormous background of photons emitted 
in jet hadronization, photons from jets are attributed to the jets and a direct photon, 
which is isolated from any jets, should be identified. Therefore, being quite different from 
J/ipiT) inclusive hadroproduction, the hadroproduction of J/ip{T) associated with a direct 
photon will be an important process to investigate the production mechanism of heavy 
quarkonium and impose new constraints on the contribution of COM to heavy quarkonium 
hadroproduction. 

We thank Gong Bin for helpful discussions. This work is supported by the National Nat- 
ural Science Foundation of China (No. 10775141) and by the Chinese Academy of Sciences 
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TABLE I: The partial cross section at LO and NLO with different cut conditions for transverse 
momentum pt but fixed rapidity cut condition \y\ < 3 for both J/ip{T) and the direct photon, 
(units: pt{GeV), (7(nb)) 



7 


Pt 


"lo 


" NLO 




T 

Pt 


T 


T 

^NLO 


Kr 


> 0.0 


> 3.0 


34 


68 


2.0 


> 3.0 


2.1 


1.5 


0.71 


> 1.5 


> 3.0 


34 


34 


1.0 


> 3.5 


1.7 


0.94 


0.55 


> 3.0 


> 4.5 


5.5 


6.0 


1.1 


> 5.0 


0.82 


0.46 


0.56 


> 5.0 


> 6.0 


1.2 


1.2 


1.0 


> 7.5 


0.24 


0.16 


0.67 



-1 

10 














NLO LHC (T) \ 




LO LHO (T) \ - 


-2 


NLO LHC (J/+) 


10 


: LO LHO (JA) \ - 




. , , . 1 . , . . 1 , , , . 1 . . . . 1 , , , . 1 . . . . 1 , , . . 1 , . . . 1 , . . . 1 , \ 



01 23456789 
P, cut of 7 (GeV) 



FIG. 6: The partial cross section of the J/ip{T) production associated with a direct photon as a 
function of pt cuts of the photon at the LHC. Here we fix the low bound of the J/'^(T)'s pt at 
lOGeV 
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